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Abstract 

The air density fluctuations in the plumes of fully- 
expanded, unheated free jets were investigated 
experimentally using a Rayleigh scattering based 
technique. The point measuring technique used a 
continuous wave laser, fiber-optic transmission and 
photon counting electronics. The radial and centerline 
profiles of time-averaged density and root-mean-square 
density fluctuation provided a comparative description 
of jet growth. To measure density fluctuation spectra a 
two-Photomultiplier tube technique was used. Cross- 
correlation between the two PMT signals significantly 
reduced electronic shot noise contribution. Turbulent 
density fluctuations occurring up to a Strouhal number 
(Sr) of 2.5 were resolved. A remarkable feature of 
density spectra, obtained from the same locations of jets 
in 0.5^ M^1.5 range, is a constant Strouhal frequency 
for peak fluctuations. A detailed survey at Mach 
numbers M = 0.95, 1 .4 and 1 .8 showed that, in general, 
distribution of various Strouhal frequency fluctuations 
remained similar for the three jets. In spite of the 
similarity in the flow fluctuation the noise characteristics 
were found to be significantly different. Spark schlieren 
photographs and near field microphone measurements 
confirmed that the eddy Mach wave radiation was 
present in Mach 1 .8 jet, and was absent in Mach 0.95 jet. 
To measure correlation between the flow and the far 
field sound pressure fluctuations, a microphone was kept 
at a distance of 50 diameters, 30° to the flow direction, 
and the laser probe volume was moved from point to 
point in the flow. The density fluctuations in the 
peripheral shear layer of Mach 1.8 jet showed significant 
correlation up to the measurement limit of Sr = 2.5, 
while for Mach 0.95 jet no correlation was measured. 
Along Ae centerline measurable correlation was found 
from die end of the potential core and at the low 
frequency range (Sr < 0.5). Usually the normalized 
correlation values increased with an increase of the jet 
Mach number. The experimental data point out eddy 
Mach waves as a strong source of sound generation in 
supersonic jets and fail to locate the primary noise 
mechanism in subsonic jets. 


I. INTRODUCTION 

The work presented in this paper addresses the issue of 
noise generation from supersonic free jets. The closing 
of the High Speed Civil Transport (HSCT) program has 
pointed out the technological advances needed in this 
arena for the success of a future supersonic, civilian 
aircraft. In addition, more stringent government and 
community requirements for aircraft noise is creating a 
further challenge to reduce noise emission from the 
current subsonic fleet. So far various experimental 
efforts depend on a "cut-and-try" approach, where one 
primarily looks for a change in the far field sound 
spectrum after making a new nozzle. A more science- 
based approach is to look for a change in the noise 
sources and therein lies the difficulty of the current 
experimental and computational techniques. In spite of 
the significant advances in the computational techniques, 
handling complex nozzle geometries and practical 
Reynolds numbers have remained a challenge. The 
theoretical background for noise emission by supersonic 
jets has been established in the past through the work of 
Lighthill (1954), Phillips (1960), Lilley (1972) and others. 
Although the recent DNS simulations (Freund 1998, 
Colonius et al. 1997) have successfully quantified the 
Lighthill stress tensor and die quadrupole sources, albeit 
for a simple geometry and very low Reynolds number, the 
same cannot be said of die experimental work. Direct 
measurement of these terms has remained an unattainable 
goal. Even at an elementary level, simple turbulence 
statistics for compressible flows are difficult to measure. 
Smits & Dussauge (1996), in a recent book on turbulent 
shear layers, write: "The form of spectrum in compressible 
turbulence is still unknown." The present work at first 
takes on die task of measuring turbulence spectra, 
specifically density fluctuation spectra, using a Rayleigh 
scattering based technique. This is followed by a direct 
correlation study "between the cause and the effect (Lee & 
Ribner, 1972)", that is, between the turbulent density 
fluctuation and far field noise. The goal is two fold: first, 
to provide a reliable and accurate database that can be used 
to validate computational aeroacoustics codes; second, to 
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locate and determine the relative strength of various sound 
sources experimentally. 

With the realization that large-scale vortical 
structures play a significant role in the development of 
free shear layers (Crow & Champagne, 1971, Tam & 
Burton 1984) a route different from the acoustic analogy 
approach (Lighthill 1954) began to emerge. 
Experimental work on low Reynolds number jets by 
McLaughlin, Morrison & Trout (1975) and Morrison & 
McLaughlin (1979) provided sufficient proof that in 
supersonic jets large organized vortical structures are 
responsible for a significant part of sound generation. 
This approach is more 'physically realizable' and has 
been successful in predicting noise emitted from 
supersonic jets in the form of 'Mach wave' emission 
(Seiner, Bhat & Ponton 1994). However, the noise field 
of subsonic jets and the high frequency part of the 
acoustic spectra could not be explained satisfactorily. 
The present work deals with both subsonic and 
supersonic jets. The experimental results were not 
strongly tied to any theory. It should be pointed out that 
the noise sources can be expressed solely in terms of 
density through the left hand side of the Lighthill's 
equation: t^p/dt 2 - ao 2 V 2 p = ^Tjj/SxjSxj (where p is air 
density, ao is sound speed, and Ty is the turbulence stress 
tensor). The properties of large vortical structures, 
required in the instability wave based approach, can also 
be expressed through density p(x,r,cp,t) = 
p"{r}exp[i{ax+n<p-cot}] (where p" is the eigenvector, 
r,<p,x are cylindrical coordinates and a, n and © are 
wavenumbers). 

The lack of experimental data in compressible 
jets is primarily due to various problems in using 
traditional experimental tools. In the past, dynamic 
measurements in compressible jets were attempted using 
hot-wire probe or Laser Doppler Velocimetry (LDV). The 
success of using hot-wire by McLaughlin, Morrison & 
Troutt (1975) and Morrison & McLaughlin (1979) is 
attributed to the very low Reynolds number operation of 
their jet. At higher Reynolds numbers hot-wires are prone 
to breakage. At mixed subsonic and supersonic operating 
conditions die analysis of hot-wire signals become 
intractable. LDV, along with a high seeding rate, has been 
used (Lau 1981, Jiang & Sislian 1998, among others) to 
determine velocity fluctuation statistics. However, various 
problems associated with seed particles following the 
flow, severe biasing problem in turbulent shear flows and 
difficulties in measuring spectra cast a doubt on such data. 
Similar problems associated with seed particles are 
encountered in Particle Image Velocimetry (PIV). In the 
present program a laser-based, non-intrusive Rayleigh 
scattering technique that depends on light scattering from 
air molecules is used. Since no seed particles are used, 
many problems associated with LDV and PIV can be 
overcome. The work of Wilson & Damkevela (1970) 


needs special attention, as they measured mean-square 
density fluctuation using a non-intrusive crossbeam 
technique. The technique involves cross correlating 
deflections from two orthogonal laser beams passing 
through the jet When the beams intersect, the cross- 
correlation provides a mean square density fluctuation at 
the point of intersection. Some of their results are 
compared with similar data obtained from the present 
study. The crossbeam technique was also used to correlate 
density fluctuations to far field sound pressure fluctuations 
(Damkevala etal. 1973, among others). For this purpose, 
the beams were separated in space and information along 
the complete path of the beams was integrated. Moreover, 
the assumption of homogeneous turbulence, used in this 
technique, is invalid in a jet All of these cast doubts on the 
validity of their results. 

Based on Proudman's (1952) analysis, the above 
mentioned casuality approach of relating flow fluctuations 
to the far field noise was taken by many researchers. In 
this method, the velocity, density or pressure fluctuations 
in die jet were correlated with the sound pressure 
fluctuation measured by a microphone. The primary 
problem of all of these efforts is with the flow fluctuation 
measurement The noise produced by intrusive hot-wire 
probe (Seiner & Reethof, 1974) or microphone placed 
inside die flow (Meecham & Hurdle, 1973, Hurdle, 
Meecham & Hodder, 1974) contributed most of the 
correlation in some experiments. Later on LDV was used 
by Schaffer (1979) and Richarz (1979) among others. 
Schaffer wrote down die vast number of approximations 
and assumptions needed to relate the experimental data to 
theoiy. Nevertheless, an important issue of accuracy in 
velocity spectra measurements using LDV (Schaffer 
seeded only the primary jet not the entrained flow) was not 
answered. Before attempting to measure correlation, it is 
natural to expect that die individual components would be 
measured accurately. This important step, however, was 
skipped in all of the earlier efforts. A significant part of the 
present work concentrates on this aspect before attempting 
a correlation measurement. It is interesting to note 
Professor Ffowcs Williams' (1973) evaluation report on a 
host of experimental papers at the 1973 AGARD 
conference on Noise Mechanism: "Many papers presented 
at die Specialist Meeting dealt with different aspect of 
source location but none of them faced up to the difficult 
issue of interpreting or speculating on the degree of 
ambiguity that must inevitably be contained in the 
experimental results." A brief review of experimental 
techniques used for noise source locations can be found in 
Grosche (1979). 

A part of the present work is an improvement 
and extension of our earlier effort to measure density 
fluctuation in supersonic jets (Panda and Seasholtz 1999b, 
henceforward referred to PS 1). A modified optical setup 
and data reduction technique, presented in this paper, is 
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found to minimize the effect of electronic shot noise in 
spectral information. Instead of measuring the scattered 
light by a single photo-multiplier tube (PMT), it is split 
into two parts and measured using two PMTs. Since the 
shot noise generated by individual detectors is 
uncorrelated, a cross-correlation between die two brought 
out die desired flow fluctuation information. The jet 
facility was also modified to suit flow-acoustics 
correlation measurements. 

n. EXPERIMENTAL SET-UP 

Experiments were performed at NASA Glenn Research 
Center using three different nozzles (one convergent and 
two convergent-divergent) operated at Mach numbers, M 
= 0.95, 1.4 and 1.8. The convergent-divergent nozzles 
were designed by the method of characteristics and their 
geometries were reported in Panda & Seasholtz (1999b). 
All nozzles were 25.4 mm in exit diameter. The jet facility 
used a continuous supply of unheated compressed air. The 
facility was located in a large test chamber, which was not 
anechoic per se, but acoustic absorbent material was 
placed around the vicinity of the nozzle and in the ceiling 
and walls of the test cell to minimize reflection. Two 1/4 
inch microphones with die 'protection grid' removed were 
used to measure sound pressure fluctuation spectra. The 
coordinate system used throughout this paper is shown in 
fig. 1. The traversing system allowed flow measurements 
in the horizontal x-y plane. The microphones were placed 
in the vertical x-z plane to minimize the effect of reflection 
from the large optical lenses and beam traps. Whenever 
possible, various parts of the optical train, air supply duct 
and the jet facility were covered by 1/4 inch thick 
polyurethane foam. However, the uncovered optical 
surfaces in the jet vicinity remained acoustic reflectors. 
The Rayleigh scattering system is somewhat elaborate and 
the following provides a description. 

A. Laser Rayleigh scattering: 

It is known that in Rayleigh scattering, the intensity of 
scattered light is proportional to the number density of gas 
molecules and in turn, is related to the local air density (p). 
For the present experiment light intensity was measured 
using photomultiplier tubes and photon counting 
electronics. Therefore, the number of photoelectrons (N) 
counted over a time interval At can be expressed as: 

N = kg At , (1) 

where k is a proportionality constant to be determined 
through calibration. Figure 2 shows a schematic of the free 
air jet facility and the Rayleigh scattering setup. The 
Rayleigh scattering based point measurement technique 
used a continuous wave (CW) laser beam focused to a 
point in the flow, with the molecular scattered light 
collected and measured using PMTs. The green laser light 
(532-nanometer wavelength) from a CW, single 
frequency, frequency-doubled Neodymium Vanadate 


(Nd:YVO<) laser was transferred to the jet vicinity by a 
0.365mm core diameter, multi-mode optical fiber and 
focused to a probe volume. Light scattered by the air 
molecules were collected at a 85° scattering angle and 
focused on the face of another 0.55 mm core diameter 
multi-mode optical fiber. The combination Of focusing and 
imaging optics makes the probe volume length of 1 .03 
mm. The complete transmitting and receiving optics were 
mounted on a x-y traverse that allowed the probe volume 
to be moved automatically over a plane. A somewhat more 
detailed description of the optical system can be found in 
Panda & Seasholtz (1999c). The scattered light, collected 
by the receiving fiber, was transmitted to a quiet 
neighboring room away from the noisy environment. Here 
the collected light was collimated and then split into two 
equal part parts by a thin-film beam splitter. Each of the 
beams was refocused into individual photomultiplier tubes. 
Photon counting electronics were then used to measure 
light intensities. The advantage of the photon counting 
approach over the conventional measurement of analog 
PMT output (Pitts & Kashiwagi 1984; Gouldin & Halthore 
1986; Graham, Grant & Jones 1974) is a clearer estimate 
of measurement uncertainty due to electronic shot noise. 
The counting was performed over a series of contiguous 
time bins of specified interval. The maximum number of 
bins that can be used at a time was 16384. Usually 
multiple sets of data were collected and passed on to a 
Personal Computer. The timer cycle for bin width was 
supplied externally through a programmable signal 
generator. All analog signals (from microphone amplifiers 
and pressure transducers indicating plenum and ambient 
conditions) were digitized and passed to the same 
computer. The data collection process was automated to 
move the laser probe volume from point to point in die 
flow field, perform die photon-counting process and 
collect the time history. When needed microphone and 
plenum transducer signals were also collected 
simultaneously. 

The success of a Rayleigh scattering system 
depends on the use of clean, particle free air, minimization 
of stray scattered light and providing a stable, vibration 
free environment for some optical components. The last 
requirement is not as critical as for velocity and 
temperature measurement (where a Fabry-Perot 
Interferometer is used). Since, the facility was built to 
measure velocity and temperature in addition to air 
density, special arrangements were made to minimize the 
effect of loud jet noise. The primary jet was supplied with 
unheated, compressed air filtered to remove all dust 
particles. A clean co-flow through a coaxial nozzle of 200- 
mm diameter surrounded the primary jet. An external air 
filter and air handling system were installed to produce 
this co-flow. The coaxial flow ensured dust-free air for the 
entrained flow. However a few particles were unavoidable. 
The second problem was stray reflected laser light. Due to 
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the special orthogonal arrangement of the optical setup, the 
primary laser beam intersected the nozzle block when the 
probe volume was close to the nozzle exit Between the 
nozzle exit and 1 .3 diameter downstream of the exit the 
stray background light overwhelmed the Rayleigh 
signature and no measurement was possible. The third 
problem was condensation during the operation of Mach 
1 .8 nozzle. The primary supply air was dried to a dew 
point of -50° Celsius and moisture condensation was 
absent in the primary air jet. However, the entrained room 
air contained a significant amount of moisture, which 
condensed in the shear layer and caused significant stray 
scattering. An effective solution was found in operating 
die Mach 1.8 jet for a couple of hours continuously before 
starting the data acquisition. The jet facility was located in 
a closed test ceil and a couple of hours of operation purged 
the moist room air and significantly reduced condensation. 

B. Density calibration 

A calibration process was necessary to determine the 
proportionality constant between the photon count rate and 
the air density, as well as to determine the residual stray 
light collected through the optics. The calibration was 
performed in the unheated plume of a convergent nozzle 
operated in the Mach number range of 0 to 0.99. At each 
operating condition the photon arrival rate was counted 
over a second duration and the jet density is calculated 
using isentropic relations. Subsequently, a straight line was 
fitted through the data to determine die proportionality 
constants a and b. 

N = (aq +b)AL (2) 

The additional constant b is needed to account for the 
room light and stray scattered laser light. Figure 3 shows 
sample calibration curves. Since two counters were used, 
two sets of calibration constants ai, bi and a 2 , b? were 
calculated. The calibration was performed over a density 
range smaller than that encountered in the jet plumes. This 
is not of concern since the fundamental linear relationship 
between die molecular number density and intensity of 
scattered light holds true at any gas density. 

C. Time average density and root-mean-square density 
fluctuation measurement: 

The instantaneous flow density p is divided into a time- 
averaged part p, and a fluctuating part 
q =~f +q’ (3) 

Photoelectron counting (Nj, i=l,2,..n) over a large number 
(n = 65536) of contiguous time bins covering more than 1 
second time duration was used to measure the two parts. 

The mean density was related to the mean count rate N^ 
by 

p~ = — ( Nav -b) , where Nav = . (4) 

a n 


All density statistics presented here had to deal with the 
electronic shot noise arriving from the photo-multiplier 
tube. (No confusion should be made between ‘jet noise’ 
and ‘shot noise’; die former is an acoustic phenomenon 
while the later is of photo-electronic in nature. Any optical 
intensity measurement is inherently affected by shot 
noise). Even when the incident light is of constant intensity 
(no density fluctuation), the rate of photoelectron emission 
by a PMT shows significant variation referred to as 
statistical photon count noise or ‘shot-noise’. This noise is 
random in nature and follows Poisson's statistics. The 
density fluctuations cause the collected light intensity to 
vary and the joint statistics can be expressed through 
Mandel's Formula (Saleh & Teich 1991). An important 
result of the Mandel's formula is that the variance of the 
photon count c^n is a sum of the variance of shot noise, 
cr 2 sh, and the variance of the light power fluctuation, crp 2 : 

<4 = ■ (5) 

The variance of photon count was calculated from the 
measured data as: 


n 

<4 = 1 




( N i-N g J i 
Nr 1 


( 6 ) 


Another important result of Poisson's statistics is that the 
mean square of shot noise is equal to the time average of 
all counts: 


< Ty/j (7) 

This allows for the light power fluctuation measurement: 
a P 2 = ct n 2 - Nay. It is straightforward to show that the 
equivalent variance in the density fluctuation, p' 2 , can be 
determined through die calibration constant a, and the time 
width At, of the individual bins: 




2 _ <?> 

n 2 


( 8 ) 


The root-mean-square density fluctuation, p^ was 
calculated as 
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(9) 


To reiterate, photon counting was performed over 65536 
time bins. The average (equation 4) and mean square 
(equation 6) were calculated; the former was subtracted 
from the later, the mean square fluctuation was determined 
through equation 8, and finally a square root of mean 
square provides rms density fluctuation. 

The fundamental source of error in the time 
averaged density data is once again from the shot noise. 
The relative uncertainty in the measurement of N is 
determined through equation 7 as 
1 


G’sh —tt ~ r 


N a 


= No 


( 10 ) 
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For the present experiment the count rate was high: 
between 8 and 12 million per second. Due to the small 
time width, the count accumulated in the individual bins 
was small, yet averaging over the large number of bins 
reduced the uncertainty to <.05%. Since, the contribution 
to uncertainty from the fundamental noise source was 
veiy low, that from a host of secondary sources became 
prominent. Occasional particles were unavoidable and 
their passage led to an increase in count. There were 
very fine oil droplets, perhaps picked up from the air 
compressor that caused a small difference in scattering 
intensity between the supply air and the cleaned ambient 
air. These secondary noise sources are difficult to 
quantify; the absolute density numbers are found to be 
repeatable within ±1% of their quoted values. 

The primary error source in rms data was due to 
a small deviation from Poisson's statistics (equation 6) in 
the PMT and photon counting electronics. Data obtained 
from quiescent clean air condition (steady light 
intensity) shows that the average count deviated slightly 
(±2%) from the mean square calculation. Since, the 
difference between the mean square and mean count was 
used to calculate the root-mean-square density 
fluctuation, the above deviation produced a positive bias 
error of about 4%. The rms fluctuation measurement is 
susceptible to unsteadiness in the laser intensity, which 
fortunately, was very small. The second major source of 
error was due to the passage of occasional particles 
through the probe volume. As the probe volume was 
moved from close to the nozzle exit to 14 diameters 
downstream, the number of particles increased 
progressively from a few per second to the order of 1 00 
per second. No measurements were performed beyond 
x/D = 14. Large particles passing through the probe 
volume were easily detected by examining the bin to bin 
variation, and removed by neglecting counts above 5 
times the rms value. This procedure, however, could not 
account for the small increases caused by particles 
passing through the vicinity of the beam waist. The bias 
error increased progressively from the nozzle exit to 
farther downstream as a larger number of particles were 
entrained into the jet. No data were taken beyond a 
downstream distance of 14D. 

D. Measurement of density fluctuation spectra 
The straightforward route to measure density fluctuation 

spectra is to take a one-sided Fourier transform of the 

sequence of photon counts N(t), 

(11) 

0 

and then subtract a constant level that corresponds to the 
average shot noise contribution: 


Shot noise floor = 

Nav (12) 

no. of frequency bins in spectrum 

This was the route taken in the earlier effort (Panda & 
Seasholtz, 1999). The method provided a crude estimate at 
the energetic low frequency part, and was unable to 
resolve the high frequency part of the spectrum. Shot-noise 
contributes to a fixed floor and a superimposed 
randomness to spectrum. The latter could not be removed 
through the subtraction process. In the present technique 
the collected light was split into two nearly equal parts and 
measured with two PMTs. The simultaneous photoelectron 
counting produced two series of dataNi(t) and N 2 (t). The 
average counts from each of the time series were 
subtracted: N,(t) = Nj(t) - N lav , N 2 (t) = N 2 (t) - N^ and a 
cross-correlation, 

00 

N\N^(t ) = J N\(t) N' 2 (t+z) dt , (13) 

—OO 

removed die shot noise contribution for the most part. A 
cross spectral density function was used 

^N\N' 1 (f) = N\N' 2 (r)e llj ft dr 


= F(N\)-F'(N' 2 ) 

to calculate the density fluctuation spectra: 


/>(/) = 


F(N' i )-F m (N' 2 ) 

a } a 2 (At) 2 


(15) 


The F indicates Fourier transform and F* complex 
conjugate of the transform. Usually a long record of either 
262,144 (for higher Mach number jets) or 524,288 (low 
Mach number jets) data was made from multiple segments 
of 16384 data strings. The latter is the maximum number 
of contiguous counts delivered by the photon counters. 

Two such records were collected from two PMTs. The 
Welch method of modified Periodograms (1967) was used 
to calculate die cross-spectral density. Each long record 
was divided into small segments of 5 12 data points. The 
adjacent segments were overlapped by 50%. The modified 
periodograms of corresponding segments from the two 
PMTs were calculated and then used to determine local 
estimates of cross-spectral density. All local estimates 
were averaged to obtain the final cross-spectral density. 

A source of experimental uncertainty in the 
spectral Hata is due to aliasing from the unresolved part of 
fluctuation spectrum. The photon counting process does 
not allow for an external anti-aliasing filter as used in the 
traditional signal analysis. However, the counting process 
sums up all fluctuations over the bin duration which acts as 
a special type of low pass filter. To illustrate, a 10 
microseconds bin width effectively yields a sampling rate 
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of 100 kHz. and the spectral information is resolved up to 
50 kHz. A model of the summing up process shows that 
the contribution from fluctuations higher than the Nyquest 
range (50 kHz.) progressively falls down to zero at the 
sampling frequency (100 kHz.). In addition, the low 
frequency, higher energy containing part of spectra is 
always resolved. The less energetic high frequency part 
does not affect measurements significantly. Therefore the 
aliasing error is expected to be small. 

E. Measurement of flow-sound cross correlation: 

A V* inch diameter microphone was kept at a fixed 
position in the far field and the laser probe volume was 
moved from point to point in the flow field to determine 
correlation between the flow density fluctuations and 
sound pressure fluctuations. The microphone was kept at 
either 30° or 90° to the flow direction (Fig. 1) and at a 
radius of 50 nozzle exit diameters. The microphone signal 
was measured through an analog-to-digital (A/D) 
converter and the Rayleigh scattered light through a 
photon counting process. An external timer clock was used 
to time synchronize the acquisition processes. Towards 
this goal a signal generator was programmed to produce 
square waves of specified time period. The A/D 
conversion occurred at each rising edge of the square 
waves, which also marked the beginning of individual 
time bins for photon counting. A time bin ended at the 
following rising edge. To verify time synchronization of 
this process, a synthetic signal was measured and 
compared. The synthetic signal was used to drive an 
electrostatic actuator and a photo-diode. The electrostatic 
actuator provided an excitation signal to the microphone 
diaphragm and the photo-diode, placed in front of the 
collection fiber, produced light intensity modulation. 
Satisfactory synchronization, observed over all frequency 
ranges tested, provided confidence in the data acquisition 
system. 

The microphone was calibrated using a piston 
phone and the calibration constants were used to convert 
the voltage signal to instantaneous sound pressure 
fluctuation p'. The cross correlation density function 
between die flow density and sound pressure fluctuations 
are calculated using the sequence of photon counts and the 
microphone outputs: 

0000 

FitpCD** 4 \\N<(t)p'(t+z) e -** dtdx. 

(16) 


» w a At 

The cross-spectral density has real and imaginary parts, 
F e y 09 = Re F p .p.(f) - rim F p . p tf) , (17) 

which are used to determine the magnitude and coherence 
function. 


00 

Fni'p'O) 


F t y (f\ = V Re2 + Im 2 



Fn&O) 

Fn\n\ 0) Fp'p'O) 


( 18 ) 


Since the collected light was split and measured with two 
counters, the cross-correlation was performed two times: 
between microphone signal and either one of the two 
series of counts. Finally, an average of the two results was 
calculated. The cross-correlation minimizes shot noise in 
the photon count data, but a noise floor persists. To lower 
the noise floor more than half a million (524,288) data 
points were collected for each data string. The Fourier 
transform used die segmenting and averaging process as 
outlined in the previous subsection. Unlike the calculations 
for density spectra, longer segments of the data string 
(4096) were used to account for die propagation time 
delay. The sound waves radiated from the flow 
fluctuations require finite time (distance of 
microphone/sound speed) to reach the far field 
microphone. Therefore, the time signature of the 'effect 1 
(sound pressure fluctuation) is expected to lag behind the 
time signature of 'cause' (turbulent fluctuations). From a 
consideration of the sampling rate and microphone 
distance it was determined that for the present case only 
first few (<5%) and last few data strings were affected by 
the time delay. By changing the length of data strings it 
was confirmed that the final correlation remained 
independent of this parameter. 


in. RESULTS AND DISCUSSION 
The choice of the three Mach number conditions (0.95, 1 .4 
and 1 .8) was deliberate. The purpose was to cover a range 
where both subsonic and eddy Mach wave radiation 
mechanisms produce jet noise. The schlieren photographs 
and tiie microphone spectra presented in Figure 4 shows 
that this goal has been attained. (In Fig. 4a a previously 
obtained schlieren photograph for Mach 0.99 jet is used.) 

It is known that in a supersonic jet, when turbulent eddies 
attain a speed faster than the ambient sound speed, a Mach 
wave emission process begins (Lowson & Ollerhead 1968; 
Bishop, Ffowcs Williams & Smith 1971; Papamoschou 
1997) The Mach waves are ballistic shock waves attached 
to the supersonic eddies. The difference in the schlieren 
photographs for the three jets are primarily due to the 
inception of the eddy Mach wave emission process. Traces 
of these waves are visible in the photograph for the Mach 
1 .4 jet, while a stronger radiation pattern is visible in 
Machl.8. There is a wide spread, between 0.89Uj to 0.6Uj, 
where Uj is the jet centerline velocity, in the convective 
velocity of turbulent eddies reported by various 
researchers. The estimated convective velocity (Table I) is 
always subsonic for Mach 0.95 jet; the Mach 1 .4 jet is at 
the borderline, while eddies in Mach 1 .8 jet are expected to 
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attain supersonic speed. To determine the frequency of 
emission a microphone was kept in the near field at the 
position marked by ® in the schlieren photographs. The 
respective microphone spectra are shown in Fig. 4(d). 
These spectra confirm the eddy Mach wave radiation for 
Machl.4 & 1.8 jets and their absence in Mach 0.95 jet. 
There is a common sharp hump at a lower Strouhal 
number (Sr = fD/Uj, f = frequency in Hertz, D nozzle exit 
diameter, Uj jet centerline velocity) around 0.2 in the near 
field spectra of fig. 4(d). The hump is associated with the 
large vortices present at the end of jet core. Zaman (1986) 
reported that the footprints of the vortices are seen along 
the entrained flow streamlines. Eddy Mach wave emission 
dominated die high frequency part of the sound spectra 
and peaked around Sr = 1 .5. The far field sound pressure 
spectra are presented in Fig. 4(e). The 30° angle was 
chosen to coincide with the peak sound radiation direction 
for cold jets as reported by Yu & Dosanjh (1972). It 
should be pointed out that the ripples superimposed on the 
microphone spectra are telltale signs of reflection caused 
by large optical components placed around the jet. 
Nevertheless, the primary difference in spectra between 
the subsonic and supersonic conditions is seen as an 
increase in die overall level and a more energetic high 
frequency part. In a nutshell. Fig. 4 demonstrates that the 
three Mach number jets used in this experiment covered 
both subsonic and Mach wave sound emission regimes. 

Before going into a discussion of density 
measurement, it is perhaps worthwhile to present some 
important aspects of viewing a jet flow through this 
parameter. In a low speed unheated air jet, the density 
variations are negligibly small. As the jet velocity is 
increased (or if heated) the density difference between the 
ambient air and the jet core increases. For the present 
unheated case increasing Mach number caused an increase 
in cooling which, in turn, increased die density of the 
primary air jet The turbulent density fluctuations are 
caused by simple mixing between the ambient and primary 
jet fluid, as well as by the inertial effects (local 
acceleration and deceleration) of flow. The former perhaps 
contributes more to the density fluctuations; although, the 
inertial effects are expected to be the source of acoustic 
radiation. Table I shows the density difference present in 
the current experimental conditions. The accuracy of the 
current technique to measure density fluctuations improves 
as the density difference increases. Therefore, 
measurements made from higher Mach number conditions 
are expected to be more accurate. The experimental data 
are non-dimensionalized by the difference between the jet 
centerline and the ambient density, (pj-p a ). The time- 
averaged data ( p) were non-dimensionalized as ( p- 
P»)/(Pj-Pa). Here pj is the jet centerline jet density 
calculated from isentropic relations and p a is the ambient 
density. The parameter is unity at core and drops to 0 as 
the ambient condition is reached. The fluctuating density 


data were also normalized by (pj-p a ), i.e, Pim/(PrPa)- A 
minor disadvantage of this normalization is that the 
uncertainty in the finally presented data increases with a 
decrease of the jet Mach number. The same uncertainty of 
±1% in the absolute measurement of time-averaged 
density manifest as an error of ±5% in the non- 
dimesionalized presentation of data from Mach 0.95 jet 
and at ±1 .5% for Mach 1 .8 jet. 

A. Time averaged data 

The centerline variation of time averaged density is shown 
in Fig. 5(a). The C-D nozzle surveys were performed at a 
pressure ratio that produced the weakest shocks in the 
core. It is to be pointed out that a completely shock free 
plume is never realized from the C-D nozzles at supersonic 
operating conditions. A literature survey of earlier data 
confirms this observation. Supersonic jets are wave-guides 
and a small manufacturing irregularity or even the shear 
layer turbulent fluctuations will produce standing waves in 
the form of weak shocks. A least shock operation point is 
determined by changing the operating pressure ratio and 
measuring the shock strength. For the present nozzles the 
least shock operation is measured to occur at M = 1 .395 
and 1.795 which are close to the design conditions of M = 
1.4 and 1.8. The residual weak waves produce significant 
density modulations that persist far beyond the potential 
core (Fig. 5a). The inset schlieren photographs confirm the 
presence of the weak shocks. Downstream from the nozzle 
exit the spatial distance between the peaks and valleys 
reduces progressively, perhaps as a result of the lowering 
of local Mach number. For the present unheated jet, 
centerline density decays to the lower ambient value. Jiang 
& Sislian (1998) used Rayleigh scattering to measure 
centerline decay of heated jets where density increased to 
the ambient value. 

The root-mean-square density fluctuations of Fig. 
5(b) shows that the turbulent fluctuations remain low over 
the first 6-10 nozzle diameters, then grow exponentially at 
the end of the potential core, and finally taper down 
slowly. As the potential core becomes longer with an 
increase of Mach number, the peak positions of the rms 
density fluctuation also moves further downstream. The 
peak fluctuations is about 0.22 (pj - p a ). Note that, similar 
to the time averaged data, the fluctuating quantities are 
also normalized by the difference between centerline and 
ambient density. The relative uncertainty level in the 
fluctuation data is higher due to the shot noise 
subtraction process described in the earlier section. In 
addition to a 4% bias error, there is a random uncertainty 
of ±4%. As mentioned earlier, the number of particles 
passing through the probe volume increased progressively 
and particularly affected data points obtained beyond the 
potential core. The measured data from this region are 
expected to be biased towards even higher values. The 
root-mean-square density' fluctuations can be nearly zero 



under two situations. First if the turbulence fluctuations are 
very low, as in the potential core of the jet, and second if 
the flow velocity decays to the incompressible regime, 
such as in the far-field of jet development. In the later 
situation, there would be a significant velocity fluctuations 
while density fluctuations would be insignificant. 

The jet spreading, due to the growth of shear 
layer, is shown in the time-averaged radial surveys of Fig. 

6. The corresponding time-averaged velocity and 
temperature data, measured using spectrally resolved 
Rayleigh scattering technique, can be found in Panda & 
Seasholtz, 1999a. The radial profiles demonstrate a 
progressive growth of the shear layer. The higher Mach 
number M = 1 .8 jet has the expected trend of a slower 
spreading rate compared to the Mach 0.95 case. The 
deviations from top-hat distribution, in the measured 
profiles from close to the nozzle exit, are due to the 
presence of the weak shock-cells. The radial profiles of 
root-mean-square density fluctuation are shown in Fig. 

7. At the closest measurement station from the nozzle 
exit, the fluctuations peak in the shear layer around y/D 
= ±0.45. The fluctuation profile at x/D = 2 shows a fully 
turbulent shear layer with a quiescent core. The peak 
locations progressively approach the centerline, and 
onward from about 10D for Mach 0.95 jet and 12D for 
Mach 1.8 jet the centerline is the location for peak 
fluctuations. There is very little information on 
experimental measurement of scalar turbulent 
fluctuations in the available literature. One exception is 
the crossed-beam measurements of Wilson and 
Damkevala (1970). The cross-correlation between 
signals from two perpendicular, intersecting beams 
provided measurements of root-mean-square density 
fluctuations. In a Mach 0.6 jet they reported peak 
pnns/(pj-pa) ~ 0.2 that agrees well with the current data. 

By presuming that the ambient pressure is 
imposed everywhere in the jet plume (i.e., local pressure p 
= pamb) it is possible to calculate time-averaged 
temperature from density. The inverse of die density_ 
profile becomes the temperature profile T = p an t/( pR), 
where R is the universal gas constant (Jiang & Sislian, 
1998). The above assumption may not be valid in a 
supersonic stream, where weak waves can produce 
pressure variation. It is also known that even in a subsonic 
jet local static pressure can deviate from the ambient 
pressure (Hussain & Clark 1977, Zaman 1986). Therefore, 
no attempt is made in the present paper to calculate local 
temperature. 

C. Density fluctuation spectra 
Figure 8 shows the outcome and improvements attainable 
using 2 PMTs over single PMT. The width of time bins 
used in the photon-counting process was 10 microseconds. 
Following die Nyquest criterion, fluctuations up to 50kHz. 
were resolved. Figure 8(a) shows spectra of the individual 


PMT signals. Following equation 12 the dashed line shows 
the shot noise contribution. In Panda & Seasholtz (1999b) 
the fluctuating part above the shot noise floor was 
presented. The improvement caused by cross-correlating 
the two PMT signals is shown in Fig. 8(b). The density 
fluctuation spectrum is better defined; the improvements 
are especially prominent at the high frequency end. 

Figures 8(c) and 8(d) show, respectively, the coherence 
and phase relation between the two PMT signals. The 
coherence function varies between 0 to 1 . A value of 0 
would be measured if both signals were solely due to the 
electronic shot noise, i.e., there were no correlation 
between the two signals. A value of unity indicates perfect 
correlation. Therefore, the measured plot of the coherence 
function shows relative importance of the density 
fluctuation to the fluctuations from the electronic shot 
noise at various frequency ranges. Since, the same density 
fluctuations were measured by both PMTs, there should be 
no phase difference between the two signals. Figure 8(d) 
confirms this expectation and increases confidence on the 
measurement technique. 

The cross-correlation process cannot eliminate 
noise from the density spectra; it lowers the noise floor. 

The finite number of data points used for the correlation 
and other electronic noise determine the final floor. An 
estimate of this noise floor is shown in Fig. 9 where the jet 
is turned offbut the clean co-flow is present Since there 
are no density fluctuations the spectral content is mostly 
due to the residual noise. The passage of occasional 
particles produces an additional contribution at the low 
frequency end (Seasholtz & Panda 1999), seen as a small 
ramp in Fig. 9. The particular data set for Fig. 9 was 
obtained 10 diameters away from the nozzle exit where 
particle entrainment rate was relatively high. Spectra 
obtained at no flow condition and close to the nozzle exit 
do not show die low frequency contribution. Nevertheless, 
all density spectral description should be compared with 
the no-flow baseline and the part above the baseline is 
truly the desired spectrum. A comparison between Figs. 9 
and 8(b) shows that at the high frequency end the 
measured fluctuations approach the residual noise floor. In 
addition to the noise base, die density spectra have a 
random uncertainty, the extent of which is visible in the 
randomness superimposed on the basic spectral shape of 
fig. 8(b). 

The measurement technique was applied to 
determine the Mach number effect on density fluctuation 
spectra. For this purpose the laser probe volume was kept 
at a fixed position in die shear layer and the plume Mach 
number was varied between 0.5^M^1.8. The resultant 
spectra are shown in Fig. 10. When frequency values were 
normalized to Strouhal number an interesting fact 
emerged: die peak frequency of turbulence fluctuations 
was always found to occur at Sr = 0.65 at the particular 
measurement station. Hot-wire measurements at low Mach 
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number jets, M < 0.5, show a similar trend. The fact that 
this trend extends to supersonic conditions is shown for the 
first time. Another observation is that the spectra become 
fuller as Mach number is increased. The reason may not be 
completely due to a general increase of turbulence level, 
but instead may lie with an increase in the signal to noise 
ratio. As the Mach number is increased the difference of 
density between the ambient and the jet flow increases. 

The higher die difference, the better the signal to noise 
ratio. This in turn lifts a larger part of the spectrum above 
the noise floor. 

Figure 1 1 presents a comparison between the 
density fluctuation spectra in Mach 0.95 and 1 .8 jets. Parts 
(a) and (c) were measured along the peripheral shear layer 
and parts (b) and (d) were along the centerline. The data 
from the shear layer shows a gradual shift in the spectral 
peak from Sr = 1 .0 to 0.2 with an increase in the 
downstream distance. The trend follows velocity 

fluctuation measurements in a low speed jet (Crow & 

Champagne, 1971) and expectations from hydrodynamic 
stability analysis. The latter establishes that an increase in 
shear layer thickness with downstream distance leads to 
increased amplification of progressively longer 
wavelength and lower frequency Kelvin-Helmholtz 
instability waves. The spectra measured along the 
centerline (Figs. 1 lb and 1 Id), however, differ from the 
velocity fluctuation measurement at low speed. Unlike the 
low speed measurements there is no definite hump at the 
passage frequency of large organized structures. All 
spectra peak at Sr = 0 for Mach 0.95 jet. For Mach 1 .8 jet, 
weak humps are visible close to the nozzle exit, yet further 
downstream the peak occurs at Sr = 0. The primary 
observation of Figure 1 1, however, is the overall similarity 
in the spectral shapes, both along the centerline and the 
shear layer, between the two Mach numbers considered. 

This similarity is further explored in fig. 12 where 
fluctuations occurring in the individual Strouhal 
frequencies are plotted from a large data set The density 
fluctuation spectra were measured at a large number of 
points ( 8 to 12 axial bylO radial points), and the mean- 
square fluctuations occurring at desired frequency bins 
were separated. For a uniform comparison, the square root 
of the mean-square fluctuations were non-dimensional tzed 
by the difference of density between the jet core and the 
ambient density. In other words, from equation 1 5 the 
plotted variable can be written as: 



Pj-Pa Pj-Pa 

where the frequency bin fo corresponds to the desired 
Strouhal number. The process is die same as applied to the 
earlier root-mean-square data. Figure 12 shows similarity 
in the distribution of turbulent fluctuations in the three 
Mach number jets considered. Fluctuations in the highest 


Strouhal number occur close to the nozzle lip and in the 
shear layer. The position of the peak fluctuations 
progressively moves downstream as the lower Sr are 
considered. The Sr = 0.25 shows highest fluctuation at the 
end of the potential core and at y/D ~ 0.3. The fluctuations 
occurring at even lower Strouhal number, e.g. Sr = 0.06, 
peaks close to the centerline and further downstream. The 
primary change visible in die fluctuation pattern is an 
overall stretching with an increase in the Mach number. 

The slower growth of the instability waves and the 
resultant slower mixing process is the reason behind this 
stretching process. Returning to the sound pressure 
fluctuation data of Fig. 3, it is seen that the noise 
characteristics of the three Mach number jets are 
significantiy different. This is in contrast with the 
similarity in the turbulent fluctuations. This led to the 
following correlation measurement between the flow 
fluctuations to die far field noise. 

D. Flow-sound correlation: 

The correlation data presented in this paper are in terms of 
normalized cross-spectral density in stead of the traditional 
cross-correlation coefficient. The primary advantage is the 
additional knowledge of frequency distribution that can be 
gained from such a representation. The intermediate steps 
towards the calculation of cross-spectral density between 
density fluctuations and sound pressure fluctuations are 
illustrated in Fig. 13. Part (a) shows the sound pressure 
fluctuation spectrum in a dimensional form (Pascal , no 
conversion to die common deciBel scale is performed). 
Part (b) shows the density spectrum measured using two 
PMTs. Part (c) presents the magnitude of the cross-spectral 
density function (equation 16). Finally, part (d) shows the 
coherence function obtained by normalizing cross-spectral 
density by the microphone spectrum and the spectrum of a 
single PMT output. The coherence function y provides a 
measure of linear dependency between the flow 
fluctuations and the sound fluctuations. All correlation 
data presented in the rest of the paper are in terms of y. 
Note that the two PMT correlation process, used for 
density fluctuation measurement, was not used for flow- 
sound correlation. Instead each of the PMT outputs was 
separately correlated with the microphone signal and an 
average between the two results was finally presented. 

Before proceeding further an estimation of noise 
floor in the calculations of coherence function is 
necessary. Similar to the density fluctuation data the cross- 
correlation data are affected by electronic shot noise. The 
shot noise is independent of the microphone signal; 
therefore, the cross-correlation process is expected to 
remove it However, a residual base level was always 
retained. To estimate this level, the jet was turned off 
(keeping the co-flow on) and the-coherence function 
between the microphone and the PMT signals was 
measured. Figure 14 shows the base noise level obtained 
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from this noflow data. A comparison between die flow 
and noflow data, respectively Figs. 13(d) and 14, shows 
that coherence function measured beyond 20kHz is 
completely due to the base noise level. However, for 
fluctuations slower than 20kHz, the measured coherence is 
significantly above the noise base. In other words, flow 
fluctuations up to 20kHz are radiating noise in the far field. 
Since, a larger coherence also implies more efficient 
radiation process, density fluctuations occurring below 
1 0kHz are more effective (y » 0.06) than fluctuations 
between 10 - 20kHz. range, where coherence 
progressively falls to the noise floor (y as 0.002). In 
addition to the noise floor there is a convergence 
uncertainty that contributes towards the randomness seen 
in the coherence function plot. 

A comparative study of correlation between 
turbulence fluctuations occurring along the jet centerline 
and sound pressure fluctuations at 0 = 30° and R = 50D for 
the 3 Mach number jets is shown in Fig. 15. The 
microphone was kept fixed while the laser probe was 
moved from point to point along the centerline. Note that 
the ordinates are in a linear scale as opposed to the 
logarithmic scale used for Figs. 13 and 14. The individual 
plots are shifted and their new bases are marked by dotted 
lines. The no flow results at the bottom of each figure 
show the noise floor. In Mach 0.95 jet (Fig. 1 5a) no 
correlation was measured till the laser probe was moved 
beyond x/D = 6. Fluctuations occurring at lower 
frequencies (up to Sr = 0.2) only are found to be 
correlated. As die jet Mach number was increased, higher 
frequency fluctuations, measured inside the potential core, 
began to radiate. For example, measurable coherence at Sr 
= 0.8 was seen for Mach 1 .8 jet (Fig. 1 5c). In addition, the 
levels were higher and a new trend in the peak correlation 
emerged. Similar to the density fluctuation spectra, the 
peak correlation progressively occurred at lower frequency 
as the laser probe was moved downstream. The peak 
coherence was around 0. 1 and measured in the region 8 < 
x/D <, 12. 

The flow sound correlation was also measured 
using a microphone at 0 = 90°. This microphone was 
significantly more affected by acoustic reflection from the 
uncovered optical elements, beam dump and the 
uncovered part of the jet facility. In turn, the flow-sound 
correlation was also affected. Figure 16 presents some data 
obtained from Mach 1 .8 jet using a 90° microphone. The 
broad modulations in data seen especially in coherence 
plots from x/D = 12 and 14 are believed to be due to 
acoustic reflection. A comparison with corresponding data 
from the 30° microphone, Fig. 1 5(c ), shows a significant 
decrease in the level and spread of correlation (note 
stretched ordinate for the 90° case). Additional data 
measured for Mach 0.95 and 1 .4 jets (not shown) showed 
an even larger decrease. In fact, no correlation was 


measured between the density fluctuations occurring 
anywhere in Mach 0.95 jet and the sound pressure 
fluctuations at the 90° location. 

The correlation data for fluctuations from the 
shear layer are presented in Fig. 17. The microphone was 
fixed at 0 = 30° and the laser probe was moved along y/D 
= 0.45. A problem with the circumferential positioning of 
the laser probe needs to be mentioned. The shear layer is 
present all around die jet periphery and the correlation 
values are expected to be dependent on the circumferential 
position of the probe volume (Richarz, 1979). Referring to 
Fig. 2, the traversing arrangement allowed flow 
measurements only in the horizontal x-y plane. Since the 
optical elements for the Rayleigh set up were also kept in 
this plane, acoustics reflections were at a maximum. To 
minimize reflections, microphones were always located in 
the vertical x-z plane. The highest correlation is expected 
when both the flow fluctuations and sound pressure 
fluctuations are measured in the same x-z plane. Therefore, 
the top part of shear layer needed to be accessed. Towards 
this end the complete nozzle facility was lowered by 
0.45D: the radial position for the probe volume. The 
comparative study of Fig. 17 shows significant difference 
between the three Mach number conditions. The coherence 
function lies in die noise floor for all measurement stations 
in the subsonic Mach 0.95 jet; indicating no sound 
radiation from density fluctuations occurring in die shear 
layer. A very small rise above the noise floor is seen in 
Mach 1.4 jet, particularly at x/D = 6 and 8, indicating very 
weak radiation. Finally, Mach 1 .8 jet shows significant rise 
above the noise floor, indicating a strong sound radiation 
that can be correlated with flow fluctuations happening at 
as high a frequency as the measurement limit: Sr = 2.5. 
When these shear layer data are compared with data 
measured along the centerline (Fig. 15), additional facts 
emerge. Only low frequency density fluctuations (Sr< 0.8) 
along the centerline radiate noise while the high 
frequencies are primarily radiated from the shear layer. 

Note that the shear layer data are presented over a larger 
frequency range than the centerline data. 

The implications of all of the above observations 
are not clear at die present time. It remains to be seen how 
various jet noise theories compare in predicting the above 
trends. This section on flow-sound correlation started with 
the observation that density fluctuation spectra in all three 
Mach number jets are similar. It is now clear that the same 
cannot be said of sound radiation to the far field, which 
explains the differences in the schlieren photographs and 
the far field acoustic spectra of Fig. 4. In Mach 0.95 jet 
only a small part of density fluctuation spectra (Sr < 0.3) 
shows correlation with sound pressure fluctuations. The 
Mach 1.8 jet, on the other hand, shows strong correlation 
up to Sr = 2.5. The eddy Mach wave emission process is 
present in Mach 1.8 case and is absent in Mach 0.95 jet. 

The Mach 1.4 jet stays somewhat in between. For the 
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Mach 1.8 jet the flow sound correlation follows the pattern 
of hydrodynamic instability waves. Close to the nozzle 
exit, the thin shear layer facilitates growth of the high 
frequency organized vortices and it is that part of density 
fluctuation which is found to radiate sound. Downstream 
from the nozzle exit the shear layer progressively thickens 
and lower frequency organized vortices become more 
energetic. The flow-sound correlation also follows suit. 
The 'eddies’ in the eddy Mach wave emission are therefore 
nothing but the instability waves. Since, the instability 
waves are well predicted by the linear stability analysis, 
the success of this approach in predicting supersonic jet 
noise is now apparent (Tam & Burton 1984). 

The subsonic Mach 0.95 jet presents a different 
story. The same hydrodynamic instability waves are 
equally present in the subsonic jet, yet little correlation is 
measured between the flow and sound fluctuations. The 
density fluctuations occurring in the lip shear layer are 
uncorrelated with either the 90° or 30° microphone. Only 
weak correlation values are measured for fluctuations 
below Sr = 0.3 from the end of the potential core. The 
present technique fails to locate the primary noise source 
for subsonic jets. At this point, it should be pointed out that 
the normalized correlation coefficient (also expressed as 
coherence) shows the linear dependence between the flow 
fluctuations and radiated noise. Source (one frequency of 
flow fluctuations radiating a different frequency sound) 
and propagation non-linearity are not accounted for. It is 
possible that in a subsonic jets sound generation is 
primarily a non-linear process. The same arguments 
perhaps can be made for sound pressure fluctuations at 90° 
to the jet axis; no correlation with the flow fluctuations are 
measured beyond Sr = 0.5 at any Mach number condition. 

IV. SUMMARY AND CONCLUSION 
This paper presents an improved technique to measure 
air density fluctuation from Rayleigh scattered light, an 
application of this technique to map turbulent fluctuation 
spectra and finally relationship between the flow density 
fluctuations and the far field noise. The flows under 
study are unheated plumes of a subsonic Mach 0.95 and 
two supersonic, Mach 1 .4 and 1 .8, jets. 

The fundamental superiority of the Rayleigh 
scattering technique above the traditional hot-wire. Laser 
Doppler Velocimetry or Particle Image Velocimetry lies 
in its non-intrusiveness and its ability to obtain 
information directly from the gas molecules (no seed 
particles are used). A major problem of using the 
Rayleigh technique is in the weak scattering process. 

The light intensity measurements are significantly 
affected by electronic shot noise. A significant part of 
the paper is devoted to present the experimental 
technique, various means to reduce contribution from 
electronic shot noise, and uncertainty estimates. The 
Rayleigh scattered light from a point in the flow was 


measured using photo-multiplier tubes and photon 
counting electronics. The counting was performed on a 
large number of contiguous bins of small time duration. 
An average of all counts provided a measure of time- 
averaged density. Subtraction of the average from the 
mean square provided measurement of mean-square 
density fluctuations. The spectral information of density 
fluctuations was gathered by dividing the collected light 
into two parts, measuring individual intensities using 
two sets of PMT and photon counting electronics and 
finally cross-correlating the two time signals. Since, 
electronic shot noise from the two PMTs is uncorrelated, 
the cross-correlation process significantly reduces its 
contribution. Finally, acoustic pressure signal from a 
microphone, placed in the far noise field, was correlated 
with the Rayleigh light to determine the portion of the 
density fluctuations radiating in the far field. A single 
PMT was sufficient to measure the Rayleigh signal as 
cross-correlation with the microphone signal reduced the 
shot-noise contribution. The noise floor, in all of the 
above measurements, was determined from no flow 
conditions (jet turned off). 

Time averaged density measurements, obtained 
from radial and centerline surveys, confirmed expected 
trends such as a slow down of jet spread with an increase 
in Mach number. Nominally, x/D =1.5 was the closest 
position to the nozzle exit where a radial survey could be 
made. The root-mean-square fluctuation, normalized by 
the difference between jet and ambient density p m /(Pj-pa). 
was found to peak in the shear layer and fall to the noise 
floor in the quiescent core. The peak fluctuation was 
measured to be about 0.22 and was in agreement with 
previous measurements of Wilson & Damkevala (1970). 
The centerline surveys showed the existence of weak 
shock-cell structures inside the pressure-matched 
supersonic jets. The time averaged schlieren photographs 
confirmed this observation. 

The density fluctuation data presented in this 
paper are believed to be the first such measurements in 
supersonic jets. The sources of fluctuation are due to 
mixing between die jet and the ambient air, and inertial 
effects. The data provides a footprint of the turbulence 
fluctuations present in the jets. In order to compile a 
database for computational fluid dynamic simulations, 
density spectra were measured at a large number of points 
in the flow. When frequency values were normalized to 
Strouhal number, and spectral density values are divided 
by the difference of jet and air density, the fluctuations 
appeared remarkably similar in all three jets. Even the 
distribution of various Strouhal fluctuation components 
was similar. The primary difference was a spatial 
stretching with an increase in the Mach number. The 
distributions followed expected behavior of Kelvin- 
Helmholtz instability waves (Morris 1977). The high 
Strouhal frequency fluctuations were mostly energetic in 
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the initial thin shear layer. As the shear layer thickened 
further downstream, the spectral peaks moved to lower 
Strouhal numbers. 

Interestingly, the noise characteristics of the three 
jets, measured using a microphone and visible in spark 
schlieren photographs, were very different. The eddy 
Mach wave radiation process was present in the supersonic 
Mach 1.8 jet while it was absent in the subsonic Mach 0.95 
jet; the Mach 1.4 case showed weak Mach waves. To 
determine the sound sources a 'causality approach' was 
followed where the density fluctuations were correlated 
with far field sound pressure fluctuations. The non- 
intrusive, particle free experimental technique provided 
more accurate and reliable measurements than any earlier 
effort. Most of the correlation measurements were 
performed with a microphone placed 50 diameters away 
and at the peak noise angle of 30° to the jet axis. The data 
were expressed in normalized cross-spectral density 
(coherence) that provided frequency information. In Mach 

I . 8 jet when the laser probe was moved axially 
downstream along the lip shear layer (y/D = 0.45) the peak 
of flow-sound correlation is found to progressively move 
from die maximum resolved frequency of Sr = 2.5 to a 
lower Sr = 0.5. The peak coherence value also increased. 

A similar trend is seen for correlation from various points 
along the centerline. All of these point out that die Kelvin- 
Helmhohz instability waves were direcdy responsible for 
eddy Mach wave radiation. The scenario changed 
significantly for the subsonic Mach 0.95 case, the density 
fluctuations in the peripheral shear layer were found to be 
uncorrelated with die far field sound. Relatively weak 
correlation, only in the 0 £ Sr ^ 0.3 range, was found from 
the end of the potential core and along the jet centerline. 
The same Kelvin-Helmholtz instability waves are present 
in the subsonic jets, yet die primary sound radiation 
mechanism is different and could not be identified in the 
present study. 
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Microphone 


Operating conditions Table I. 


Sp. heat ratio, y = 1.4; Total temp., To - 300°K; Ambient density, p a 
1.16 Kg/m 3 , Ambient sound speed, a = 347 m/s 


Nozzle type 

Convergent 

M 1.4 CD 

M 1.8 CD 

Minimum shock operation at 
M = 


1.395 

1.795 

Jet Velocity Uj (m/s) 

316 

411 

486 

Estimated eddy convection 
speed U c (m/s) = 0.6Uj - 0.89Uj 

190-282 

247-366 

292-433 ] 

Jet density pj (kg/m 3 ) 

1.36 

1.6 

1.89 

Difference (pj -pj (kg/m 3 ) 

0.2 

0.44 

0.73 

Frequency for Sr = 1 (kHz) 

12.4 

16.2 

19.1 



Fig 1. The coordinate system 
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Photon arrival rat© (per sec) 


Anechoic box 



7x106 



Fig 3. Calibration curve of Rayleigh scattered light to measure air density. Straight lines show least square fit. 
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Fig 4. (a), (b), (c) Spark schlieren photograph of fully expanded jets at indicated Mach number conditions; (d) pressure 
fluc tuati on measured at x/D=4, y/D=2 (shown by 0 in schlieren photos); (e) sound pressure fluctuations at 50D and 
30° to the flow direction. 
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Distance from nozzle exit, X/D 


Fig 5. Centerline variations of (a) time-averaged density and (b) root-mean-square density fluctuations for indicated 
Mach no jets. The insets (c) and (d) are time averaged schlieren photographs. 
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0.0 0.5 1.0 1.5 

Mean density, _ElEl . 

(Pj-Pa) 

Fig 6. (a), (b), (c) Radial profiles of time-averaged density measured at indicated axial stations and Mach no 



0.0 0.1 0.2 0.3 0.4 

rms density, . Pr ~ . 

(Pj-Pa) 

Fig 7. (a), (b), (c) Radial profiles of root-mean-square density fluctuations measured at indicated axial stations and 
Mach no conditions. 
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_2 

0 1 2 3 4 5x1 0 4 

Frequency in Hz 

Fig 8. Cross-correlation between two PMTs to obtain 
density fluctuations spectrum. Data from Mach 1.8 jet 
at y/D = 0.48 and x/D=2.4. (a) Spectra of individual 
PMT signals (b) cross-spectrum between the two 
providing density spectrum, (c) normalized cross- 
spectrum (d) cross-spectral phase. 


0 0.5 1.0 1.5 2.0 2.5 3.0 

Strouhal no, Sr 

Fig 10. Density fluctuation spectra from a fixed position 
(x/D=2.5, y/D=0.45) in the jet but at different Mach 
numbers. The vertical dashed line locates spectral peaks. 
Individual spectra are separated by a multiplication factor 
of 2. 
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Density fluctuation, (kg/m3)2 


X/D = 7.0 


(c) M = 0.95 


Lij ‘ ! J 


X/D = 12 


(d) M = 0.95 

s-« 


_ „ „ . „ Stpouhal no, Sr Sbouhalno,Sr 

»g 11. Density fluctuation spectra measured at indicated axial stations and along (a) /b) v/D=0 48 and /hi /hi 
centime of Mach 0.95 and 1 .8 j*. Individual speed* are separated by a muldpli^ fcL of 2 ® 
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(c) M = 1 .8 
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Distance from nozzle, x/D 

Fie 12. The distribution of density fluctuations (northed by pj - P.) 
different Mach numbers. Contours are at intervals of 0.002. 


at the indicated Strouhal frequencies for three 





Frequency, Hz 


Fig 13. Cross-correlation between flow density 
fluctuations and sound pressure fluctuations in Mach 1.8 
jet; laser probe at centerline and x/D=10, microphone at 
far field x/D = 50, 0 = 30°. (a) Sound pressure spectrum, 
(b) density spectrum, (c) cross spectrum and (d) 
normalized cross spectrum (coherence). 



Fig 14. Normalized cross-spectrum between density and 
sound fluctuations at no flow condition (jet turned off). 
Laser and microphone positions are as in Fig. 13. 
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Sr Sr Sr 


Fig 15. Normalized cross- spectrum between density and sound fluctuations from 3 different jets. The microphone was 
fixed at x/D = 50, 0 = 30° and the laser position was moved at indicated axial positions along centerline. 



Sr 


Fig 16. Normalized cross-spectrum between density and sound fluctuations in Mach 1.8 jet. The microphone was fixed 
at x/D = 50 , 0 = 90 ° and the laser position was moved at indicated axial positions along centerline. 
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Fig 17. Normalized cross-spectrum between density and sound fluctuations from 3 different jets. The microphone was 
fixed at x/D = 50, 0 = 30° and the laser position was moved at indicated axial positions along shear layer, y/D - 0.45. 
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